arXiv: 1509.02523v2 [astro-ph.CO] 1 Dec 2015 


Draft version December 2, 2015 

Preprint typeset using D-T^X style emulateapj v. 01/23/15 


LARGE OPACITY VARIATIONS IN THE HIGH-REDSHIFT LYa FOREST: 

THE SIGNATURE OF RELIC TEMPERATURE FLUCTUATIONS FROM PATCHY REIONIZATION 

Anson D’AloisioV Matthew McQuinn 1 , & Hy Trac 2 
Draft version December 2, 2015 

ABSTRACT 

Recent observations of the Lya forest show large-scale spatial variations in the intergalactic Lya 
opacity that grow rapidly with redshift at z > 5, far in excess of expectations from empirically 
motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in 
the ionizing background, but found that this required either extremely rare sources or problematically 
low values for the mean free path of ionizing photons. Here we report that much - or potentially 
all - of the observed excess likely arises from residual spatial variations in temperature that are 
an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity 
fluctuations generated in this way depends on the timing and duration of reionization. If the entire 
excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending 
but extended reionization process that was roughly half complete by z ~ 9 and that ended at z ~ 6. 

In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the 
most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. 
This correspondence potentially opens a new observational window into patchy reionization. 

Subject headings: dark ages, reionization, first stars — intergalactic medium — quasars: absorption 
lines 


1. INTRODUCTION 

When the first galaxies emerged ~ 100 — 500 million 
years after the Big Bang, their starlight reionized and 
heated the intergalactic hydrogen that had existed since 
cosmological recombination. Much is currently unknown 
about this process, including what spatial structure it 
had, when it started and completed, and even which 
sources drove it. The Lya forest provides one of the 
only robust constraints on this process, showing that it 
was at least largely comp lete by z ~ 6, when the Universe 
was one billion years old dFan et a l. 2006: Gallerani et~ahl 
l2008t iMcGreer et, all 1201511 . This Letter argues that 
there exists another, potentially groundbreaking signa¬ 
ture of reionization in the Lya forest data. 

The amount of absorption in the Lya forest can be 
quantified by the effective optical depth, r e g = — In (F)l, 
where F = exp[— T^ ya (x)} is the transmitted fraction of a 
quasar’s flux, (...)l indicates an average over a segment 
of the forest of length L, and TL ya (a:) is the optical depth 
in Lya at location x along a sightline. The optical depth, 
TL yQ (a;), scales approximately as the H I number density, 
which after reionization scales as T _0 7 A^/r. Here, T 
is temperature, Ab(x) is the gas density in units of the 
cosmic mean, and T(a;) is the H I photoionization rate, 
which scales with the amplitude of the local ionizing ra¬ 
diation background. 

Observations of high -2 quasars show a steep increase 
in the dispersi on of r B ff among coeval forest segm ents 
around z = 6 (|Fan et, alll200fit [Becker et alJl2015j h In 
the limit of a uniform ionizing background, the well- 
understood fluctuations in Af, fall well short of producing 
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th e observed disp ersion at z > 5.5, as shown recently by 
iBecker et al.l ( 2015 ) (hereafter B2015). Previous stud¬ 
ies have attempted to explain this excess with spatial 
fluctuations in the ionizing background. The properties 
of spatial fluctuations in the background depend on the 
number density of sources and the mean free path of pho¬ 
tons, A®( f p. While A®(f 2 is well constrained at 2 < 5.2, 
being too large to yield significan t backgro und fluctua¬ 
tions for standard source models (IWorseck et al.; 201 1 i. 
B2015 showed that the excess r e ff dispersion at 2 = 5.6 
could be matched in a model where A^ decreases by a 
factor of w 5 between 2 = 5.2 and 2 = 5.6 - a time scale 
of just 100 million years. However, such rapid evolu¬ 
tion in A®® f p is inconsistent with extrapolations based on 
measurements at low er redshifts (|Becker fc Boltoi]l201.'jl 
IWorseck et al.l[2014f) . and would imply that the emissiv- 
ity of ionizing sources, in turn, increases by an unnat¬ 
ural factor of ps 5 over the same cosmologically short 
time interval. Because of these issues, B2015 specu¬ 
lated that the excess dispersion was evidence for large 
spatial variations in the mean free patl0. Alternatively, 
fluctuations in the ionizing background could have been 
enhanced if the sources of ionizing photons were rarer 
than the observed population of galaxies. However, cur¬ 
rent models require half of the background to arise from 
bright sources with an extremely low space density of 
~ 10 -6 Mpc~ 3 ([Chardin et alll2015h. This scenario is a 
possibility of current debate (e.g. lMadau &; Haardtll2015l 
but see D’Aloisio et al. in prep.). 

In this Letter, we explore a source of dispersion in 
T e ff that has so far been neglected and that, unlike 
ionizing background fluctuations, has straightforward 
implications for the reionization process itself. In 

3 See also I Davies fc Furlanettol 11201511 . which appeared after we 
submitted this paper. 
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addition to A b (x) and r(x), the Lya opacity de¬ 
pends on T(x), mainly because the amount of neutral 
hydrogen after reionization is proportional to the 
recombination rate, which scales as T -0,7 . Previous 
attempts to model Lya opacity fluctuations had not 
included the residual temperature fluctuations that 
must have been present if reionization were patchy 
and temporally extended. As ionization fronts propa¬ 
gated supersonically through the IGM, the gas behind 
them was heated to tens of thousands of degrees 
Kelvin by photoionizations of H I and He I. After 
reionization, the gas cooled mainly through adiabatic 
expansion and through inverse Compton scattering 
with cosmic microwave backg ro und (CMB) photons 
dMiralda-Escude fc Reesl 11994 [Hui fc Gnedin! 119971: 
Hui fc Haimanl 120031 : iMcQuinn fc Unton Sanderbeckl 
201 fill . Since different regions in the IGM were reionized 
at different times, these heating and cooling processes 
imprinted an inhomogeneous distribution of inter- 
galactic temperat ur es that pe r sisted after reionization 
dTrac et all 120081: iCen et all 120091 : iFurlanetto fc Ohl 
120091 : iLidz fc Mallovl 1201 Jb We will show that these 
residual temperature variations likely account for much 
of the observed dispersion in r e s at z > 5.5, and may 
even account for all of it - a scenario that would yield 
new information on the timing, duration, and patchiness 
of reionization. 

The remainder of this Letter is organized as follows. In 
f]2]and tJHl we describe our simulations and methodology. 
In we present our main results. In fj5] we offer con¬ 
cluding remarks. We use comoving units for distances 
and physical units for number densities. 

2. NUMERICAL METHODS 

To model the impact of relic temperature fluctuations 
from reionization on the distribution of r e ff, we ran a 
suite of 20 cosmological hydrodynami cs s imulati ons us - 
ing a modified version of the code of lTrac fc Penl (2004). 
The simulations were initialized at 2 = 300 from a com¬ 
mon cosmological initial density field. We used a m at¬ 
ter p ower spectrum generated by CAMB dLewis et al.l 
l2000f ) assuming a flat ACDM model with = 0.3051, 
n b = 0.04823, h = 0.68, a 8 = 0.8203, n s = 0.9667, 
an d Yh p = 0.2453 , consisten t with recent measurements 
(jPlanck Collaboration et aDl2015 ). Our production runs 
use a cubical box with side length Tbox = 12.5h -1 Mpc, 
with N ( jm = 1024 3 dark matter particles and N gas = 
1024 3 gas cells. 

In each simulation, reionization was modeled in a sim¬ 
plistic manner by instantaneously ionizing and heating 
the gas to a temperature T re i on at a redshift of zj“{£ n - 
Subsequently, ionization was maintained with a homo¬ 
geneous background with spectral index a = 0.5, con- 
sistent with recent post- reionization background models 
(|Haardt fc Madaull2012ll . Utilizing the periodic bound¬ 
ary conditions of our simulations, we trace skewers of 
length L = 50ft -1 Mpc (following the convention of the 
B2015 T e ff measurements) at random angles through all 
of the hydro simulation snapshots. Each skewer is di¬ 
vided into N x = 4096 equally spaced velocity bins of 
size Answer = aL/N x (where a is the cosmological scale 
factor), and Lya o ptical depths are computed using the 
method of iTheuns et al.l (|1998f) . Although reionization 


occurs instantaneously at within each simulation 

box, we piece together skewer segments from simulations 
with different z” e U n to model the effect of an inhomoge¬ 
neous reionization process, as we describe further in the 
next section. 

The post-reionization temperatures in the simula¬ 
tions are relatively insensitive to the spectrum of the 
ionizing background, but they are sensitive to the 
amount of heating that is assumed to occur at the 
time a gas parcel is reionized. Previ ous c alcula¬ 
tions dMiralda-Escude & Reesl 11994 iTrac et all 120081 : 
iMc.Qninnl l2012t ! have bracketed the range of possible 
reionization temperatures to T re i on « 20,000 — 30,000 
K. (We note that previous large-scale reionization simu¬ 
lations do not accurately capture Treion, as they do not 
resolve the ~ 0.3 physical kpc ionization fronts.) Thus, 
we have run two sets of ten simulations - one set with 
T re ion = 20, 000 K and the other with T re ion = 30,000 K 
- where each set contains instantaneous reionization red- 
shifts of z^ n = {6,6.5,7, 7.5,8,8.5,9,10,11,12}. This 
redshift range spans the likely duration of reionization. 
Simulations with z™f* n >12 are driven to a common 
temperature by z < 6, so they are well approximated by 
the 2(efon = 12 simulation. 

Figure]!] shows the post-reionization thermal and asso¬ 
ciated hydrodynamic relaxation of intergalactic gas (and 
its effect on the Lya forest opacity) in our simulations. 
The top-left panel shows the volume-weighted average 
gas temperature in the ten T re ion = 30,000 K simula¬ 
tions. The bottom-left panel shows this same average, 
but limited to gas cells with densities of A b < 0.3, the 
deepest voids that dominate transmission in the highly- 
saturated z > 5 forest. Even at z ~ 5, the A b < 0.3 
gas temperatures differ by up to a factor of five between 
the simulations that were reionized at different times. 
The right panels show the H I number densities (top) 
and the transmission (bottom) at z = 5.8 for the same 
skewer through our = 12 (blue/solid) and = 6 
(red/dashed) simulations. The transmission is nearly 
zero in the z)"4 n = 12 case, owing to colder tempera¬ 
ture and hence enhanced H I densities, whereas there is 
significant transmission in the z)”( 3 n = 6 case. 

3. CONSTRUCTING THE r eff DISTRIBUTION: A TOY 
MODEL 

Reionization is a process that is inhomogeneous and 
temporally extended, unlike in our individual hydro sim¬ 
ulations. Modeling the thermal imprint of patchy reion¬ 
ization on the Lya forest thus requires an additional in¬ 
gredient: a model for the redsliifts at which points along 
our skewers are reionized. In this section, we present a 
simplified toy model to illustrate how we piece together 
skewer segments from our hydro simulations, and to pro¬ 
vide insight into how the timing, duration and morphol¬ 
ogy of reionization affect the amplitude of Lya opacity 
fluctuations. 

For illustrative purposes, let us assume that the Lya 
forest is made up of segments of equal length, l c , where 
each segment is reionized at a single redshift. Let us 
further assume that the reionization redshift of each 
segment is drawn from a uniform probability distribu¬ 
tion over the interval [z s tart, £ e nd], resulting in a global 
reionization history in which the mean ionized fraction, 
xmi(z), is linear in redshift (since xhii(z) is the cumula- 
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Fig. 1.— Illustration of the thermal and hydrodynamic relaxation of intergalactic gas after reionization and of how this relaxation 
affects the Lya forest opacity. Left panels: Evolution of the volume-weighted average gas temperature in our set of ten hydrodynamical 
simulations in which reionization is instantaneous at . Different curves correspond to different . In each simulation, the IGM 

is heated to a temperature of T re ; on = 30,000 K at redshift V V n . Right panels: Two example 10/i -1 Mpc skewer segments at z = 5.8 
from the simulations with zffjV = 12 (blue/solid curve) and = 6 (red/dashed). The two segments are drawn from the same location 

in the simulation box so they have nearly identical underlying dark matter fields. The top panel shows the physical H I number density, 
while the bottom panel shows the fraction of transmitted flux in the forest. 


tive probability distribution of the reionization redshift). 

Assuming that the reionization redshift field is in the 
Hubble flow, each segment of length l c spans A N x = 
4096 x Z c /(50/i -1 Mpc) velocity bins of our hydro simu¬ 
lation skewers. For the first segment, with reionization 
redshift z±, we take the initial AN X spectrum velocity 
bins of a randomly drawn skewer that has z™f* n closest 
to Z\. For the next segment, with reionization redshift 
Z 2 , we take the next AN X velocity bins from the same 
skewer through the simulation that has z™fo n closest to 
Z 2 (note that all simulations were started from the same 
initial density field). We repeat this procedure until an 
entire 50/i _1 Mpc sightline is filled. 

In what follows, we compare these toy models against 
the B2015 measurements of the cumulative probability 
distribution function of T e ff, P(< T e g). The measure¬ 
ments considered here span z = 5.1 — 5.9 in bins of width 
Az = 0.2. We construct P(< r e ff) from 4000 randomly 
drawn sightlines at the central redshift of each bin. For 
each redshift, we rescale the nominal post-reionization 
photoionization rate of our simulations (T = 10 -13 s _1 ) 
by a constant factor, such that our model P(< r e ff) is 
equal to the observed distribution at either P(< r e ff) = 
0.15 or P(< T e ff) =0.3, depending on which value pro¬ 
vides the better visual fit. We have performed extensive 
numerical convergence tests using simulations of varying 
resolution and box size. We found excellent convergence 
of P(< T e ff) for our production runs in both box size and 
resolution (especially for our patchy reionization mod¬ 
els). 

Figure [2] shows the P(< r e ff) of these toy models for 
a range of (z s tart ; -Zend, lc), compared against the B2015 
measurements in the two highest redshift bins (black his¬ 
tograms). The curves with shorter reionization dura- 



Fig. 2.— Comparison of cumulative probability distribution func¬ 
tions of r e ff in our toy model (colored curves) compared to the 
B2015 measurements (black histograms) and to an instantaneous 
reionization model with 2*efon — 8.5 (black/dotted curves). Top 
row: Varying the coherence length, l c , assuming a uniform -^eion 
probability between z = 5.9 — 13. Bottom row: Varying the 
redshift interval over which the values are drawn, assuming 

l c = 2bh ~ 1 Mpc. 
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tions, or with smaller l c , fall closer to the black/dotted 
curves, which assume that reionization occurs instanta¬ 
neously (at = 8.5, although it does not depend 

on this choice). Figure [2] affords three insights: (1) The 
larger the coherence length, l c , over which gas shares a 
similar reionization redshift, the larger the spread in r e ff; 
(2) The width of the observed P{< T e g) can be fully 
accounted for only if reionization ended at z < 7 and 
was well underway by z ~ 9; (3) The maximal width 
is achieved by a late-ending (z e nd ~ 6) and extended 
reionization model in which large contiguous segments 
(> 12.5 h~ x Mpc) of the Lya forest were reionized at 
z > 9. In the next section, we apply these insights to 
construct more realistic models of reionization that re¬ 
produce the observed width of P(< r e s). 

4. THE EFFECT OF TEMPERATURE FLUCTUATIONS ON 
THE HIGH-REDSHIFT LYa FOREST 

4.1. Results 


We generate physically-motivated reionization red- 
shift fields using simulations bas ed on the excursion-set 


model of reionization (ESMR) (Furlanetto et al. 

2004; 

Zahn et all 20071: Alvarez et all 20091: iSantos et all 

2010: 

Mesineer et al. 201 If), which has been shown to reoro- 


duce the ionizat ion structure found in full radiative trans¬ 
fer si mulations (IZahn et al.ll2005l |20 1 it iMaiumdar et a. 1.1 
I2014T) . In particular, our ESMR uses a realization of the 
linear cosmological density field and a top hat in Fourier 
space filter to generate a realization of reionization in 
a cubical box with side 400/i” 1 Mpc, samp l ed w ith 512 3 
cells. For details of the algorithm, see IZahn et all (120071) . 
We use the simplest formulation of the ESMR, with two 
free parameters: (1) The minimum mass of halos that 
host galactic sources of ionizing photons, M m i n ; (2) The 
ionizing efficiency of the sources, (. We tune ( to ob¬ 
tain a reionization history that is approximately linear 
in redshift, similar to those in our toy model. 

The ESMR calculation yields a cube of reionization 
redshifts. We construct mock absorption spectra by first 
tracing 50 h~ x Mpc skewers through this cube. We then 
piece together spectrum segments from our suite of hydro 
simulations, much like in our toy model, except here we 
match to the reionization redshifts along the ESMR 
skewers. This process does not account for how 
correlates with density, an effect that will underestimate 
(making our models conservative) the width of P(< T e g), 
which we address shortly. 

The left panel of Figure [3] shows a 2D slice through 
our fiducial reionization redshift field in which reioniza¬ 
tion spans z = 6 — 13. In the right panels, the blue/solicl 
curves show the fraction of transmitted flux along four 
different sightlines through this field, selected to span 
a range of r e g. The red/daslied curves show the corre¬ 
sponding reionization redshifts along the sightlines. On 
average, regions that reionize at later times yield more 
transmission, while regions that reionize at z > 9 re¬ 
sult in dark gaps in the Lya forest. The variation among 
these mock spectra is similar to the well-known varia tion 
seen in observations of the z ~ 6 Lya forest dFan et all 
120061) . In our interpretation, this variation reflects differ¬ 
ences in the reionization redshifts - and hence tempera¬ 
ture^ - between segments of the Lya forest. 

4 We find that pressure smoothing plays an insignificant role in 


The leftmost four panels of Figure 0] show P(< r e g) 
in three ESMR models that take T re ion = 30,000 K 
and M m ; n = 2 x 10 9 M 0 , with reionization histories 
shown in the inset of the 5.5 < z < 5.7 panel. The 
CMB electron scattering optical depths in these models 
are r es = 0.054, 0.068 and 0.080, within uncertainties of 
t he latest Planck measurement o f r es = 0.066 ± 0.016 
dPlanck Collaboration et all 120151) . These models are 
compared against the B2015 measurements (black his¬ 
tograms) and against the homogeneous reionization ref¬ 
erence model with = 8.5 (dotted curves). The 

blue/long-dashed and green/short-dashed curves corre¬ 
spond to scenarios in which reionization spans z = 6 — 10 
and z = 7—14, respectively. While these two mod¬ 
els produce significantly more width in P(< r e ff) than 
the homogeneous reionization model, they fall short of 
producing the full range of r e g required to match the 
observations. However, the fiducial z = 6 — 13 reioniza¬ 
tion model (red/solid curves) generally provides a good 
match to the measurements. A significant success of the 
fiducial model is that the observed redshift evolution of 
the P(< r e ff) width is reproduced without any additional 
tuning of parameters. (We have checked that this suc¬ 
cess also holds over z = 4 — 5, lower redshifts than those 
shown where the effect is reduced.) Indeed, there are 
not any parameters that can be tuned in our model to 
change the post-reionization evolution of the width. The 
pink shaded regions indicate the 90% confidence levels 
of our fiducial model, estimated from bootstrap realiza¬ 
tions, showing consistency with all the data aside from a 
single high-opacity point at z = 5.4 and at z = 5.6. 

The discrepancy at the highest opacities may arise be¬ 
cause our method of constructing mock absorption spec¬ 
tra does not capture correlations between temperature 
and density that should be present, since denser regions 
are more likely reionized earlier. Such correlations would 
act to increase the width of P(< r e a), as the denser re¬ 
gions around galaxies are ionized earlier in our models. 
One might naively think these correlations are small, 
because correlations between the density on the much 
larger scales of the H n bubbles and the smaller scales of 
voids in the forest should be weak, but calcul ations show 
that they may not be negligible (iFurlanetto fc Qhll2009l : 
iTrac et aLll20()8l iBattaglia et al.ll2f)13l) 

4.2. Effect of varying reionization model parameters 

The right panels of Figure [4] show the effect of vary¬ 
ing T r ei on and M min . For T reio n = 20,000 K, the dis¬ 
tribution of T e ff is somewhat narrower than the case 
with T re i on = 30, 000 K. Hotter temperatures are likely 
achieved towards the end of reionization, when H n bub¬ 
bles are larger and propagate at quicker speeds. A model 
with Treion ~ 30,000 K near the end of reionization, 
and smaller temperatures earlier on, would likely pro¬ 
duce more width than a model with T re ion ~ 30, 000 K 
at all times. However, any conclusions about T re ion 
prior to modeling the density/reionization-redshift cor¬ 
relations are premature. 

The blue/long-dashed and green/short-dashed curves 
in the right panels show the effect of varying M m j n . For 
the former atomic cooling halos (ACHs) curve, M m j n is 
set to the minimum mass required to achieve a halo virial 

generating r e ff fluctuations. 
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Fig. 3.- Left panel: Two-dimensional slice through the reionization redshift field of our fiducial excursion-set model. Right panels: 
Each panel shows a 50 1 Mpc skewer through the fiducial model, with the four skewers selected to span a range of r e ff values. The 
red/dashed curves show the reionization redshift values along the skewers, and the blue/solid curves show the Lya forest transmission 
at z = 5.8. The transmission curves are constructed by piecing together mock spectrum segments from our suite of hydrodynamical 
simulations, matching z^efon to the reionization redshifts from the excursion-set model. 




Fig. 4. Leftmost four panels: Cumulative probability distribution functions of r e ff measured in B2015 (black histograms) compared 
against an instantaneous reionization model (dotted curves) and against three excursion-set models of patchy reionization (other curves). 
The corresponding global ionization histories of these models are shown in the sub-panel. The pink shaded regions indicate the 90% 
confidence levels of our fiducial model (red curves) in the four redshift bins shown, estimated from bootstrap realizations. Rightmost two 
panels: The effect of varying the reionization temperature, T re [ on , and the minimum mass of halos that host galactic sources, M m i n , on 
P(< r e ff) at the two highest redshifts reported in B2015. 
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temperature of 10,000 K. For both curves, we tune ( to 
match the reionization history of our fiducial model (red 
curve in the inset). We find that the effects of varying 
Af m in are minor. 


5. CONCLUSION 

We have shown that residual temperature inhomo¬ 
geneities from a patchy and extended reionization process 
likely account for much of the opacity fluctuations in the 
z > 5 Lyct forest. Inhomogeneities in the ionizing back¬ 
ground may also contribute at a significant level, though 
current models in this vein have required very small mean 
free paths or extremely rare sources. We showed that 
residual temperature fluctuations alone could account for 
the entire spread of observed r e ff. A significant success 
of this interpretation is that it is able to reproduce the 
rapid growth of T e g fluctuations with redshift, despite 
having very little freedom in its post-reionization evolu¬ 
tion. In this scenario, the observations favor a late but 


Alvarez, M. A., Busha, M., Abel, T., & Wechsler, R. H. 2009, 
ApJ, 703, L167 

Battaglia, N., Trac, H., Cen, R., & Loeb, A. 2013, ApJ, 776, 81 
Becker, G. D., & Bolton, J. S. 2013, MNRAS, 436, 1023 
Becker, G. D., Bolton, J. S., Madau, P., et al. 2015, MNRAS, 447, 
3402 

Cen, R., McDonald, P., Trac, H., & Loeb, A. 2009, ApJ, 706, L164 
Chardin, J., Haehnelt, M. G., Aubert, D., & Puchwein, E. 2015, 
ArXiv e-prints, arXiv:1505.01853 
Davies, F. B., & Furlanetto, S. R. 2015, ArXiv e-prints, 
arXiv: 1509.07131 

Fan, X., Strauss, M. A., Becker, R. H., et al. 2006, AJ, 132, 117 
Furlanetto, S. R., & Oh, S. P. 2009, ApJ, 701, 94 
Furlanetto, S. R., Zaldarriaga, M., Sz Hernquist, L. 2004, ApJ, 
613, 1 

Gallerani, S., Ferrara, A., Fan, X., & Choudhury, T. R. 2008, 
MNRAS, 386, 359 

Haardt, F., & Madau, P. 2012, ApJ, 746, 125 
Hui, L., & Gnedin, N. Y. 1997, MNRAS, 292, 27 
Hui, L., & Haiman, Z. 2003, ApJ, 596, 9 

Lewis, A., Challinor, A., & Lasenby, A. 2000, Astrophys. J., 538, 
473 

Lidz, A., & Malloy, M. 2014, ApJ, 788, 175 
Madau, P., & Haardt, F. 2015, ArXiv e-prints, arXiv:1507.07678 
Majumdar, S., Mellema, G., Datta, K. K., et al. 2014, MNRAS, 
443, 2843 


extended reionization process that is roughly half com¬ 
plete by z ~ 9 and that ends at z ~ 6. 

Unlike ionizing background fluctuations, which do not 
necessarily signal the end of reionization, temperature 
fluctuations directly probe the timing, duration, and 
patchiness of this process. If most of the opacity 
variations owe to temperature, it would mean that, on 
average, the darkest > 10 Mpc segments of the 2 > 5 
Lya forest were reionized earliest, and the brightest 
segments last - a potentially powerful probe of cosmo¬ 
logical reionization. 


The authors acknowledge support from NSF grant 
AST1312724. H.T. also acknowledges support from 
NASA grant ATP-NNX14AB57G. Computations were 
performed with NSF XSEDE allocation TG-AST140087. 
We thank Paul La Plante, Jonathan Pober, Phoebe Up¬ 
ton Sanderbeck, George Becker, Adam Lidz, and Fred 
Davies for helpful discussions. 


McGreer, I. D., Mesinger, A., & D’Odorico, V. 2015, MNRAS, 
447, 499 

McQuinn, M. 2012, MNRAS, 426, 1349 

McQuinn, M., & Upton Sanderbeck, P. 2015, ArXiv e-prints, 
arXiv: 1505.07875 

Mesinger, A., Furlanetto, S., & Cen, R. 2011, MNRAS, 411, 955 
Miralda-Escude, J., & Rees, M. J. 1994, MNRAS, 266, 343 
Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2015, 
ArXiv e-prints, arXiv:1502.01589 
Santos, M. G., Ferramacho, L., Silva, M. B., Amblard, A., & 
Cooray, A. 2010, MNRAS, 406, 2421 
Theuns, T., Leonard, A., Efstathiou, G., Pearce, F. R., & 
Thomas, P. A. 1998, MNRAS, 301, 478 
Trac, H., Cen, R., & Loeb, A. 2008, ApJ, 689, L81 
Trac, H., Sz Pen, U.-L. 2004, New A, 9, 443 
Worseck, G., Prochaska, J. X., O’Meara, J. M., et al. 2014, 
MNRAS, 445, 1745 

Zahn, O., Lidz, A., McQuinn, M., et al. 2007, ApJ, 654, 12 
Zahn, O., Mesinger, A., McQuinn, M., et al. 2011, MNRAS, 414, 
727 

Zahn, O., Zaldarriaga, M., Hernquist, L., & McQuinn, M. 2005, 
ApJ, 630, 657 


